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ABSTRACT 


The electron photon transport code ITS has many applications in the physics and 
medical industries. The code uas originally interned for use in determining particle 
transport in thick materials. The code breaks down for very thin targets because the 
multiple scattering appro.ximation used to determine the electron deflection angles for 
thin steps is inadequate. A method of correction has been developed by Tom Jordan 
and Joseph Mack which combines a small angle approximation theory to the multiple 
scattering and an explicit large angle treatment based on a Poisson distribution. This 
method has been validated against several experiments with great success. The multiple 
scattering theory of.Moliere has also been incorporated into a correction scheme and 
shows good agreement with experimental data. 
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I. INTRODUCTION 


Particle transport dates back to the days of Lord Rutherford when the scattering of 
electrons through thin foils changed our view of the nature of matter. Today particle 
transport through matter is itself an industry. From determining the radiation damage 
effects on satellites in space to calculating the energy spectrum of medical x-ray ma¬ 
chines, particle scattering is an important part of the science industry. In the past only 
rigorous experimention provided scattering information for the scientist. Today with the 
aid of high speed computers and a better understanding of the various particle inter¬ 
actions, computational results can be easily obtained. Computer programs which model 
the transport of particles through matter were developed around 1968. .M. Berger and 
S. Seltzer at the National Bureau of Standards developed the first general electron and 
photon transport code called ETR,\N. From 1970 to 1981 eleven codes based on the 
ETR.AN model were developed, and in 1984 eight of these codes were combined into 
single code package called the Integrated Tiger Series (ITS). The eight codes differ in 
dimensional geometry' and two of the codes include transport in macroscopic electric and 
magnetic fields (Ref 1: p. 6J. 

Electron and photon transport computer codes such as ETR,^N were originally de¬ 
veloped to study the manner in which radiation scatters through thick materials. 
Quantities such as energy deposition and angular deflection are calculated using various 
.\Ionte Carlo schemes. Such schemes are based on condensed case histories, that is. a 
particle's energy, direction, and position are calculated at discrete intervals rather than 
continuously as the particle traverses the medium. Probability distribution functions 
based on this interval are then used to determine the state of the particle after each stop. 
This scheme has been very successful in the study of thick target particle transport. 
Recently, these transport codes have been applied to very thin foils lor transition radi¬ 
ation research. Discrepancies in the angular distribution of the transmitted electrons 
resulted which indicate an apparent break-down in the code calculation of the multiple 
scattering. 

The multiple scattering distribution is based on a substep size which is calculated 
from the particles radiation length and the target material density. Berger [Ref. 2; p. 143] 
has pointed out three advantages for small step sizes, two of which affect the angular 
distribution of the particles. If the step size is small, the majority of the scattering is 
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done within the material and boundary’ eOects need only be calculated in the partial 
substep at the escaping edge. The scattering in the partial substep will have a small 
impact on the total angular deflection so that crude approximations to the multiple 
scattering can be used for tliis region. The net angular deflection within a substep is so 
minimal that multiple scattering theories with this restriction are applicable. For thick 
materials where the number of collisions is large, these assumptions are valid. The ma¬ 
jority of the scattering is done within the boundary' of the material and the overall dis¬ 
tribution is calculated in this region. The number of collisions in the final substep is so 
small that their contribution to the overall angular distribution is minimal. However, for 
very thin materials, where the thickness is less than the substep size, the entire scattering 
distribution is determined from the crude approximation. 

The number of collisions presents another problem in the calculation of the angle 
deflection. Small angle approximate multiple scattering theories ignore tiie large angle 
calculation because if the number of collisions is large the majority of the scattering m ill 
be small angle scattering. For very thin materials the number of collisions can be so 
small that an occasional large angle scatter can make a large contribution to the overall 
angle distribution. It is therefore necessary to develop the large angle scattering dis¬ 
tribution along with the small angle profile. 

The Integrated Tiger Series of coupled photon and electron transport codes <11 Si 
were originally developed for thick material transport. When the ITS code for cylindn- 
cal geometry (CYLTRAN) was used for very thin free standing foils, the angular dis¬ 
tribution of the transmitted electrons had broader Gaussian forms then that of the 
available experimental data. The need for a correction to the multiple scattering dis¬ 
tribution for very thin materials existed, and several methods of solution were consid¬ 
ered. 

The best method of correction would be an explicit treatment of the scattering based 
on the screened Rutherford cross section whereby each electron collision deflection angle 
is calculated as the particle is stepped through the material. Although this method 
would give very accurate results the time inefliciency makes it impractical for computer 
programing. Any other treatment would require a theory of multiple scattering to de¬ 
termine the angular distribution. The five principle works on the subject are by Williams 
(Ref. 3). Goudsniit-Saunderson (Ref. 4], Moliere (Ref. 5). Synder-Scott (Ref. 6), and 
Lewis (Ref. 7), The theories at a glance seem as diverse as the random processes 
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themselves: howet er, m the limit of small angles they are essentially the same. A brief 
histoiy of the evolution of these theories will serve to amplify this point. 

W’hen charged particles are incident on a slab, the distribution of the scauered par¬ 
ticles is described by the well-known Bolizntan integro-diflerential equation. Bothe (Ref. 
8 : p. 11] showed that in the limit of small angles, this equation transforms into a 
Fokker-Planck type dillcrential equation. Although his own theory of multiple scatter¬ 
ing was flawed with inexact boundar>' conditions and appro.vimations, the Bothe- 
Fokker-Planck general form for the distribution would become the standard form for all 
small angle approximation theories which would follow. 

Williams was the first to utili 2 e the Fokker-Planck equation successfully. Altiiough 
his expression for the angular deflection showed agreement with the experimental data 
of the time, his theory would be overshadowed by the historic Goudsmit-Saunderson 
theory which was published one year later. Ihe Goudsnut-Saunderson theory was an 
exact treatment based on a Legendre potynominai expansion and was \alid for all angles. 
1 he theory was exact except for the assumption of equal path lengths. 

Almost a decade followed before another theory of multrplc scattcrii'.g was pub¬ 
lished. and that was tire theory of .\Ioherc. Mohere used small .mgle appro.'.nn.itions ta 
transfomr the standard tr.insport equation into a form which resembled dilliisicn. Ik 
then expanded the resulting expression in terms of an iteration constant to evaluate ilie 
distnbutron. The Mohere theory showed great agreement with experimental d.ita. cs- 
pecially for thin nuitcriais. and is the standard ilicorciical comparison for a!) scattentig 
expenments today. .Another small angle approximation siiortly followed by Syndcr- 
Scott which cfrcctisely derived the .Mohere integral equation from a standard dillusion 
equation. .A year later Lewis showed that the iniegro-dinerential equation lor sm.i1I 
angles conlormcd to the Synder-Scou expression. Finally Nigam. Savinderson. and ^'a- 
You Wu (Ref. 9; p. IW’I showed the equivalence of the Mohere and tlic Goudsmit- 
Saunderson theory in the limit of small angles. 

In the present paper, we arc concerned with scattering in very thin materials, w here 
the scattering angles arc small. Therefore the small angle approximation theories of 
Williams and Mohere arc appropriate. Exen with the best possible multiple scattering 
theory , the question of large angle deflections still remains. With very thin materials the 
number of collisions is so small that the occasional large angle scatter makes a signif¬ 
icant contribution and cannot be ignored. 



A method of solution has been developed by Thomas Jordan and Joseph Mack 
which combines the small angle theoiy of Williams and an explicit treatment of the large 
angle scattering. Chapter two will outline the theories of Williams, Goudsmit- 
Saunderson. and Moliere. Chapters three and four will describe the multiple scattering 
algorithm in the existing ITS code CYLTRAN and the correction scheme by Jordan- 
Mack.. 





11. THEORY 


A. WILLIAMS 

The particle scattering theorj- by E.J. Williams was the first theor>- to properly use 
the Fokker-Planck equation. As pointed out by Bothe [Ref. 8: p. 164] a particular sol¬ 
ution to the Fokker-Planck equation for multiple scattering has a Gaussian form. 
Williams states the same relation with reference to the general theory of errors. The 
difference between Bothe and Williams is that Williams limited the scattering angle to a 
finite value. In addition Williams set the cosine of the scattering angle equal to unity, 
whereas Bothe kept the cosine term in his solution of the scattcriiig cross section. 

To deterntine tlie multiple scattering distribution function Williams defined a limit¬ 
ing angle c-), such that on the average the particle would dcllcct once through an angle 
greater than while traversing the material. This can be written by setting the number 
of collisions in the back region equal to unity: 


I I (3.1) 

where P{<t)} is the coili'ion cross section. Since virtually all the scattering is dene at an 
angle less then <^, an approximation to the total distribution can be found by considering 
the deflections in this region. Williams showed that from the general theory of errors, 
the probability of scattering into an angle . of a particle due to collisions which are less 
than <t>, can be represented by a Gaussian of the form: 

(3.:i 

where a.] is the arithemetic mean value of a, and is given by: 




•^0 


d4> 


(3.3) 
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This general form of the theor}’ is a good approximation to the multiple scattering: 
however, when \N'illiams derived it he used the unscreened Rutherford cross section with 
the small angle approximations Isin'.'id ~ d , and cos ‘.'16 ~ 1 which resulted in: 


/’(</>) = 




(3.4) 


where N is the number of scattering atoms, t is the thickness of the scatterer. Z is the 
charge of the particle, M is the mass of the particle. )S = (v/c). and ^ = - ji -In 

addition, he further modified af by an approximation to take into account the screening 
of the atomic electrons. To graph his distribution Williams defined a unit of angle which 
would become the “natural" angle for several papers on multiple scattering. The angle 
is given by: 


0 




(3.5) 


Originally. Williams was concerned with thick targets, fast particles, and unit charges. 
In his second paper he modified a; to avoid second-order approximations he used in the 
original derivation. This modification e.\tcnded his theor> to thinner targets (U.(»l cmi. 

B. GOIDSMIT-SALNDERSON 

The multiple scattering theory of Goudsmit-Saunderson is often vjuotcd as an ex¬ 
act" theory and in the development of the most general form of the thccry no assump¬ 
tions or approximations are used. However, in the derivation of the collision cross 
section an assumption of small angles and equal path lengths is asserted such that the 
resulting distribution is a small angle approximation. The theory is founded in the basic 
Legendre polynomial property that the average value of any polynomial after n cvems 
is equal to the average value of the poiynonuals after one event to the n" power: 

< cos 01 >„, - < cos P.) >a, (3.01 


where 0 is the final scattered angle after several collisions, and fl, is the scattered angle 
after one collision. The total average of any Legendre polynomial can be written as: 


G, - V n (w) < P,{ cos 0,) >gy 

fImC 


(3.7) 
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The scattering distribution is then given by summing all these averages per unit solid 
angle; 


yW“-:^X^2/+l)G,/>^cos(^) (3.S) 

The collision probability can be represented by a Poisson distribution wliich has func¬ 
tional dependence on the total collision cross section p , and the thickness of the 
scatterer t: 


ll\H)~-^^,v~np\\i (3.9) 

nl 

This expression would be exact if the true path length were equal to the foil thickness, 
and for small dellcciions this assumption is quite valid. 

The evaluation of the Goudsmit-Saunderson theoiy has been developed extensively 
by Spenser (Ref. luj and Berger |Ref. 2; p. 20?J. The Goudsmit-Saunderson multiple 
scattering distribution is given by Berger as 


f ';)e\p 


(■C 


- I Gfs'ilh' Pj, c 

J 


cos t'M 


The expansion coellicients G arc gi\cn by: 


(3.Ml 





a(0)l I — P/{ cos sin 6 dO 


(3.111 


where oiO) is the screened Rutheribrd cross section which has the Ibrm; 

aid }' 




I - cos e -f 2ii )‘ 


(3.I2I 


(.3.13 


The constant xi is Moiiere's screening angle and is given by equation (3.21). To solxc 
equation (3.10). recursion relations are developed so that a large number of expansion 
coelTicients can be calculated. Berger IRcf. 2: p. 213) shows that a convenient form of 
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the distribution for random sampling can be found by replacing the Legendre 
polynomials P{ cos tu) by the relation: 

cos wj = J Pl^x) dx (3.14) 

•'cos CO 

with the recursion relations: 

Wo=l-costx> (3.15) 

//| » 'ii(l — cos'eu) (3.16) 

(/+!)///=( 2/-1) costa//,2 )//,_2 /^2 (3.17) 

This method of solution to the GouJsmit-Saundcrson theory is used in the 11S code 
s>stem for the multiple scattering within the material. 

C. MOLIERE 

.Molicrc takes a different approach in that he explicitly starts his theor>‘ with the 
assumption that all scattering angles arc small, such that the standard transport 
equation reduces to an expression which resembles ditfusion in a plane: 


1 

a 






t)<»l/ll/\ 


t3.lSi 


where <»(/) is the differential scattering cross section. N is the number cT scattering atoms 
per cni. 0, •• is a vector representing the direction of the electron before the las', 

scatter. d\ « / «// . and/?(/./) is the scattering angular distribution. .\ l ouiier 

(BesselI transformation which is gixen in great detail in Bcthc s paper produces the an¬ 
gular distribution function of the form: 


/(/,/)■[ If «/»//p(i/<») X exp[-.V/[ o{}f)/di\\-J^Ui/)}'\ t3.l'5) 

*'0 •'o 


where JXtiO) and J.Uii) itre bcsscl functions of the zeroth order. Molicrc takes the theory 
a step further by transforming his equation into an expression which depends on two 
angles 7 ’ and y]. 
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The transformed Moliere multiple scattering distribution function then takes the fol¬ 
lowing form: 


J\6)0 dO = / d/. f y dyjg{/.y) exp [ ' '*y^{ -b -1- In ! (3.2U<7l 

Jn 


(3.20i») 




y * y.cn 


(3.20c) 


O.lOd) 


The first angle is the screening angle. which describes the scattering atom and is 
derived from the specific form of the given scattering law: 

, r (//i.i'T ir i-r r 7 

" L“ JL w J'J ' 

In Moliere s original paper, he used the Tliomas-rcrmi potential for the single scattering 
law. Since this form does not contain the Born appro.ximation only elastic collisions 
against the Coulomb field of the nucleus was considered. Fano I Ref. 11: p. 1 IT] showed 
that the correction due to the inelastic collision with atontic electrons consisted of re¬ 
placing 2^ by Z(Z-^ 11 and adding an additional term to Molicre's constant b: 


b - ln( 1 * + (Z !)■'{ In [O.I6UZ" -’ (I * 3.33Zc‘//n)] - r,} 


(3.2;) 


where «• is an integral oscr an incoherent scattering function whose value is about -5.0. 
Berger (Ref. 2: p. 207). extrapolates this constant for scseral materials which are shown 
in Figure I on page lu. 

The second angle is a unit probability angle, /f. which is a measure of the foil 
thickness and states that the total probability of a single scatter at an angle greater than 
fr, is equal to one. 
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rigiire I. 


C’f vs Z < \loniic NiimlH'rK I .ino's corretiion const;int. r, verses 
utotnic miinbcr. Vuluo luivc been e\irapo)jteJ Irom Heiccr s book, 
(Ref. 2: p. 2nT|. 
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This "natural" unit of angle was first determined by Williams in his theory of multiple 
scattering. The unit probability angle is given as: 


2 4,T*V/e''Z(Z + 1 )r 

>.(= — 


(pv)- 


(3.231 


To solve the transformed equation Moliere then defines an iteration constant B„ and 
expands the integral in terms of a power series in I B„: 

B„ = b + lnB„, (3.24) 

jm de = di{f{i) + B'YiO + + ... ] (3.25) 


where 


(3.2(.') 


.Moliere s thcor\ is valid for small angles less than thirty degrees |Ref 12: p. 1256]. His 
theoiy cuii be extended to larger angles by the introduction of a multiplication factor. 
^ sin 0/6 , into the overall angular function as suggested by Bethe IRef 12: p. 1263). 
The first order approximation to the Moliere angular distribution with the correction 
just mentioned has the following form: 


, sin A .— f- 


(3.27) 


This represents the Gaussian term for multiple scattering. Higher order terms are 
sonrewhat more diflicvilt to calculate and impractical to program. Bethe indicates that 
when the scattering is Ic's than two degrees, the second term/', represents a ten percent 
correction to the Gaussian first term. (Ref. 12: p. 1260] 

Since we are concerned with only verx' thin distances where the scattering is mostly 
in the forward direction the first term will represent a good approximation to the multi¬ 
ple scattering. 
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III. ITS (CVLTRAN) 


The ITS code CYLTRAN is the transport code for cylindrical geometry. A com¬ 
plete description of the code system is provided in Appendix B. CYLTRAN uses two 
methods for determining the multiple scattering distribution of the electrons. When the 
scattering is done within the material, and within a full substep interval, a subroutine 
named ".ML’LT" is called to calculate the scattered angle of the electron. Subroutine 
MULT uses the Goudsmit-Saunderson theory to determine the scattering distribution. 
For scattering in the final partial substep at the escaping edge of the material, a sub¬ 
routine named "ANGLE" is called. ANGLE uses a crude approximation to the multiple 
scattering which is described in detail in Appendix C. Both methods hinge around the 
substep size which depends on the scattering material. The substep size is calculated by 
the equation: 


n _ Drasce 

L>subuep 


(4.1) 


where is the mean free path of the incident electron (g cm-). ISUB is the number 

of substeps per interval, and p is the material density. The scattered angle of the incident 
electron is then calculated after each substep from the Goudsmit-Saunderson distrib¬ 
ution. In a cross section generating program XGEN, the average cosines based on the 
same substep size given in the equation above, are calculated and stored in an output 
file. CLYTRAN reads in these stored cosine averages. In a subroutine named .MULT 
the angle distribution is formed using up to 240 Legendre polynominals with the cosine 
averages. Scattered angles are then drawn from the distribution for each incident 
electron after each substep. This process works fine until the final substep at the mate¬ 
rial boundary where the distance from the last full substep to the escaping edge is less 
than the substep size. In this region the cosine averages calculated by XGEN are no 
longer valid and the Goudsmit-Saunderson distribution can not be used. To determine 
the scattering angle distribution a crude approximation to the multiple scattering is used. 
The approximation is made in a subroutine called ANGLE. ANGLE is called whenever 
the distance from the last full substep to the escaping edge is less than the substep size. 
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The crude approximation is based on Williams' Gaussian expression: 


( 0 ^ \ r n -COS0) "1 

J[e)oc expf j ^ exp 1^- 5 ^ J (4.2j 

where a is calculated from the relation: 

a = (l -<cosg>) ~^ - - - (4.3) 

^substep 

The average cosines are given by the program XGEX, and SHD is the distance from the 
last full substep to the escaping edge. For thick foils this appro.ximation is reasonable 
because it has ver\’ little effect on the total distribution which is calculated from several 
Goudsmit-Saunderson calculations within the material. For very thin foils where the 
thickness is less then the step size, the entire angular distribution is determined by this 
crude approximation. It is therefore necessary to develop a better appro.ximation to the 
multiple scattering in the final partial substep. 
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IV. THE JORDAN-MACK CORRECTION 


A correction to the multiple scattering distribution which uses the Williams small 
angle approximation theory and an explicit treatment of large angle scattering has been 
developed by Thomas Jordan and Joseph Mack. The major diflerence between the 
Jordan and Williams solution is that Jordan uses the Goudsmit-Saunderson expression 
for the screened Rutherford cross section whereas Williams used the original Rutherford 
cross section. The Goudsmit-Saunderson expression is given by: 


277e^Z^ _ sin $ d9 _ 

p\^ (\-cQ^e+ 


(5.1) 


where 6 ] is Moliere's screening angle y ]. The mean square angle for single scattering is 
then defined by integrating the average angle over the cross section in the forward re¬ 
gion; 


r'irr 

0? = AV da 
*'0 


(5.2) 


where 6 „ is defined as the maximum angle such that the probability of a scatter greater 
than is exactly one. Zerby and Keller (Ref 13: p. 202) state the same relation in their 
discussion of Moliere's theorx', where they assert that the scattering angle must be re¬ 
stricted to exclude large angles for which .Moliere's theory breaks down. To accomplisli 
this, they suggest choosing an angle such that on the average one large scatter will 
occur in an angle larger than . This method was first derived by W’illiams almost 27 
years earlier. The small angle approximation is then assumed in the form: 

e^ = 2(l-cos^) (5.3) 

and the resulting mean square angle is given by: 



cos 6 ) da 


(5.4) 
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Equation (5.4) can easily be integrated into the following form: 

= 1] (5.5) 

where y\ — )il{l — + »/), = cos d„ , and »/ = Vzd\ = Vzxl ■ See Appendix B for the 

complete proof The term (1 — /i + »/) can be found from the condition that one scatter 
will occur in the back region: 


r da = I 


(5.6) 


and upon integration it yields: 


(1 = 


1 



(5.7) 


With tlie mean square angle completely determined, the Gaussian form of the distrib¬ 
ution is known; 

face CO- (5.Si 


To extend this relation to larger angles the multiplication factor sin 6/0 as suggested 
by Bethe [Ref 12; p. 126.1) for .Moliere s theory can also be incorporated into the previ¬ 
ous equation so that the final distribution function has the form: 


F{$)oc 


d 


[5.9} 


For a given electron scatter through a material distance t. the small angle defiection can 
he determined from the distribution function F{d) : however, the large angle deflection 
must also be determined. To accomplish this, a Poisson distribution function is con¬ 
structed whereby each term represents the probability of a large angle scatter; 





(5.1U) 


where m is the probability of a single event. Since the expected number of collisions in 
the back region, (6 > 6J has been set to unity, the probability of a single event is equal 
to one, (i.e. m = 1). A random number can then be drawn and the probability of each 
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individual large angle scatter subtracted from it until the remaining probability is less 
than or equal to zero. For each large angle scatter an explicit random angle determi¬ 
nation is done based on the collision cross section. The Poisson distribution can only 
be used when the scattering probability is small and constant as in the case of ver> thin 
foils. If the foil is too thick or the energy of the incident electron is too small, then the 
probability of large angle scattering will increase and Poisson statistics will not properly 
described the distribution. 

Thus the multiple scattering has been decoupled into small and large angle scatter¬ 
ing. This technique hinges around the angle parameter 6 „ , the maximum angle lor 
which one large scatter will occur at an angle greater than In the case of ultra thin 
foils or very high energy electrons it may not be possible to have a large angle scatter. 
In this case cos turns out to be greater than one! Therefore it is necessary to calculate 
the angle dellection by the explicit Poisson treatment for all collisions within a step si/c. 
For this case the probability of a single event m is equal to the total collision cross sec¬ 
tion integrated over all angles; 






V. RESULTS 


The Molicrc and Jordan-Mack correction schemes were programed into the llS 
code CYLTR--\N by inserting u call statement in subroutine ANGLE. This statement 
calls a new subroutine ANGDET (ANGle DETcrmine) which contains the correction 
schemes. When the distance from the last full substep to the escaping edge is less than 
the substep size, ANGLE is called by the main program and ANGDET j.s called by 
ANGLE. The fortran source code for both corrections is given in Appendi.x C. along 
with the existing ANGLE subroutine code. I hc input parameters to C\’LTRAN for 
each run where similar to the example input file given in Appendix B. The number of 
histories for each run was IiKt.iKKi. The statistics for each run were calculated based on 
ten batches. 

The problem of validating and comparing the existing C^ LTILAN code with the two 
correction schemes was dilTiculi due to the lack of sulTicient experimental data for scat¬ 
tering through \cr> thin loiis. Originally onh th.c historic results of Hanson [Ref. Id: 
p. o.'idi were considered. In this experiment Hanson used a lS.t>o mg enf 
(‘J.t’7 X :/.!•» \m) gold foil with 15.77 MeV electrons, figure 2 on page IS shows the 
CLYTRVN solution along with the Hanson data. Tiic substep si/e for MeV 
electrons in gold is .'2.uS mg cm- x 10 'em) so that the Hanson foil is smaller then 
the intersal subsiep. Therefore, the entire distribution is calculated in the subroutine 
.ANGLE by the crude approximation, from figure 2 on page IS it is shown that this 
approximation giscs a much wider Gaussian width then the experimental dat.i. 
figure .1 on page 19 and figure 4 on page 2<» show the corresponding .lordan-Mack and 
Molicrc methods ol correction with the Harsson data. Both correction methods show 
an overall improvement in the angular distribution. Jt is of no conscijucncc that the two 
solutions give similar results because they both consist of small angle approximate 
Gaussian forms. It will be shown later that there are regions where these methods of 
correction give dificrcni results. 

The next step was to compare these two correction methods with experimental data 
of different energies and materials. Only two other experiments for very thin foils could 
be found, and they were the papers of Kagcyama and Mozicy. (Ref. 15: p. .^4S and. 16: 
p. 647] The Kagcyama experiment was for 1.66 MeV electrons through both copper and 
lead foils, for this experiment, the low speed of the electrons ctVcciivcly increased the 


17 




Figure 2. CYLTRAN Suiutioii nilii Hanson's data for Gold: Angular Spectrum 
of transmitted electrons for 15.77 .McV incident electrons through a 
18.66 nig cm’ gold foil. 
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Figure 3. Jordan-Mack Solution with Hanson's data for Gold; Angular Spec¬ 
trum of transmitted electrons for 15.77 McV incident electrons through 
a IS.66 mg cm* gold foil. 
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Figure 4. Moliere's Solution nith Hanson's data for Gold: Angular Spectrum 
of transmitted electrons for 15.77 MeV incident electrons through a 
1S.66 mg cm* gold foil. 
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number of collisions and a broader Gaussian distribution resulted. The foil thicknesses 
wore 6.3 mg cm- (7.05 x 10 "‘cm) and 6.2 mg cm’ (5.d0 x 10“^cm) respectfully. I he cor¬ 
responding substep sbes were 14.67 mg cm’ (1.64 x 10-’em) and 8.96 rng cm’ 
(7.86 X 10-^cm). Therefore, the copper foil is thicker then the substep sbe and is seg¬ 
mented into three substeps. In the first two full subsieps, the multiple scattering dis¬ 
tribution is calculated by the Goudsmit-Saunderson theory. The third substep is less 
than the substep size so that the deflection angles are calculated by the approximation 
in the subroutine ANGLE. On the other hand, the thickness of the lead foil is less then 
the substep size so that only the .ANGLE subroutine is used for the scattering distrib¬ 
ution. Figure 5 on page 22 shows the CYLTR.AN solution for the copper case, and 
again the distribution is much broader than the experimental data. The Jordan-Mack 
method (Figure 6 on page 23), docs not show exact agreement as in the Hanson case; 
however, it does show an overall correction. For this case, the Moliere model 
(Figure 7 on page 24) showed the closest agreement to the data. 

The lead foil solutions. (Figure 8 on page 25 thru Figure 10 on page 27) show the 
same trend. 1 he CVLTR.-\N‘ solution is always broader then the experimental data and 
for low energies, the Moliere method is better than the Jordan-Mack technique. 

On the other end of the energy spectrum, the .Mozley experiment was at a much 
higher energy. 600 McV’. and consisted of scattering through very thin aluminum foils 
(2.44 mg cm’, 9.0 X lO-N'm). For this case the substep size is equal to 3.99 gem’ 
(1.4S X 10‘‘cni). so that the foil thickness is much less than the the substep size. .-\t such 
high energies (600 .MeV), the scattering is in the forward direction and not much greater 
than 0.02 degrees. This makes it very difficult to resolve the scattering and tally the 
electrons into such small angle bins. Therefore the curves appear crude. Even with 
double precision, only a bin width of 0.025 degrees was allowed before floating point 
underflow errors resulted. In spite of these computational problems, the general trends 
in the solutions can still be seen. In Figure 11 on page 28 the CLVTR.AN solution is 
given for this case, and the poor agreement to experimental data is clearly shown. The 
Jordan-.Mack solution (Figure 12 on page 29) shows great agreement with the data and 
is superior to the .Moliere case (Figure 13 on page 30). 

With the JM solution verified over a range of energies and materials, the research 
effort appeared to be winding down, until something unexpected resulted! \\'hen the half 
width values for several gold foil thickness calculations were plotted for both the 
CYLTR.AN and the JM solutions, a periodic "saw-tooth" pattern resulted (Figure 14 
on page 32). Only when the periodicity was determined to be the substep size did this 
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C^'LTRAN Solution with Kageyama's data for Copper: Angular 
Spectrum of transmitted electrons for 1.66 MeV incident electrons 
through a 6.3 mg cm’ copper foil. 
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Figure 6. Jordan-Mack's Solution with Kageyania's data for Copper: Angular 
Spectrum of transmitted electrons for 1.66 .MeV incident electrons 
through a 6.3 mg cm^ copper foil. 
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Figure 7. Moliere's Solution witli Kageyama's data for Coppen Angular Spec¬ 
trum of transnuttcd electrons lor 1.06 MeV incident electrons throuch a 
6.3 mg cm^ copper foil. 
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Figure 8. CVLTRAN Sulutioii with Kageyama's data for Lead; Angular Spec¬ 
trum of transmitted electrons for 1.6G MeV incident electrons through a 
6.2 mg cm- lead foil. 
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Figure 9. Jordan-Mack's Solution with Kageyama's data for Lead: Angular 
Spectrum of transmitted electrons for 1.66 .\leV incident electrons 
through a 6.2 mg cm’ lead foil. 
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Figure 10. Moliere's Sulution uitli Kageyama's data lor Lead: Angular Spec¬ 
trum of transmitieJ electrons for 1.66 McV .nciJcnt electrons throuch 
a 6.2 mg cm*’ lead foil. 
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Figure II. CVLTRAN Solution uith Mozley's data fur Aluminum: Angular 
Spectrum of transmitted electrons’ for 6C0 MeV incident electrons 
through a 2.44 mg cm- aluminum foil. 
















Figure 12. Jordan-Mack'.s Soliitiun «Hh Mozkv's data for Aluminum: Vn^-ular 
Spectrum of transmitted electrons for OUO MeV incident electrons 
through a 2.44 mg cm*' aluminum foil. 
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Figure 13. Muliere's Solution niili Mozley's data lor Aluminum: Angular Spec¬ 
trum of transmitted electrons for WK> MeV incident electrons through 
a 2.4d mg cm*' aluminum foil. 
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pattern make sense. In the original ll'S code, the angular distribution within a full 
subslep is calculated by the Goudsmit-Saunderson theory. If the distance between the 
last full substep and the escaping edge is less than the substep length, the subroutine 
ANGLE is called to calculated the distribution for the region using the crude approxi¬ 
mation. As shown in the theory section, this approximation is dependent on the dis¬ 
tance between the last substep boundaiy and the escaping edge. The larger this distance 
is, the larger the error in the approximation will be. However, if the foil thickness is 
exactly equal to the substep si/e, or a multiple of this size, the entire distribution will 
be calculated by the Goudsmit-Saunderson theory and the subroutine ANGLE will 
never be called. The cfTcct is such that the error from the approximation grows until the 
distance between the last substep and the edge is equal to a substep and then the sol¬ 
ution "corrects ' itself by a Goudsmit-Saunderson calculation. With this explanation, a 
better understanding of the subroutine ANCiLE is given. For the existing ITS material 
interface approximation there is a strong dependence on the foil thickness. On the otlter 
hand, the .lordan-Mack method shows ver> little dependence on this distance parameter 
SI ID. For very thick foils the error of approximation at the escaping edge is a small 
fraction of the total di'tribution which has already been established by the multiple 
substep GS calculations within the material. A semi-log plot shows this graphically 
(Figure 1.^ on page 

.An interesting result of the Jordan-.Mack solution is the apparent evaluation of 
Molicrc s iteration /?.. Upon comparison of the two fenns of the Gaussian distribution, 
(equations and (.''.do. a Jordan-.Mack constant can be formed and is given by; 

yj- - 1 

Comparisons of li, and /?. are shown in Figure 16 on page .'4 and Figure 17 on pace 

In Figure 16 the energy is held constant at 15.77 McV. while the foil thickness is 
varied from 10 microns to one meter. I his figure shows the region where the Jordan- 
Mack solution and the Moliere theory give similar results. At around eight centimeters 
B. becomes constant whereas still increases. Since the thickness of the final partial 
substep is usually much less than a centimeter this region is of little concern. Figure 17 
on page 35 however docs contain an important result. In Figure 17 on page 35. the foil 
thickness is held constant and the incident electron energy is varied from one electron 
volt to several .VIeV. Again there is a region where the two theories arc similar. At 
around 0.5 McV the two iterative constants level off. but for energies less than this, the 
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Figure 14. Hair Widths (CVLTRAN with Jordan-Mack; Linear): 1’iic input 
parameters were for 15.77 Mev incident electrons through gold foils of 


various thickness. 
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Figure 17. Iterative Constant Comparison (Moliere v«ith Jordaii-Macki: The 
input parameters were for a 18.66 mg cm* gold foil with incident 
electrons of varving energy. 
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two curves diverge. Recall that this iterative constant B is a measure of the Gaussian 
width of the scattered distribution. Since B is in the denominator of the exponential, 
when B increases the width of the distribution increases. Physically, when the incident 
electron energies arc low the number of collisions should increase and the Gaussian 
width of the distribution should increase. Therefore the value of B should increase as 
the energy decreases. The Jordan-.Mack solution appears to break down at energies less 
than 0.5 .MeV. To investigate this problem a "stand-alone" fortran program named 
JORDAN was utilised. JORD.-\N is simply a bin tally program which uses the 
Jordan-Mack correction scheme to form the cumulative distribution function from 
which the incident electron angle deflections are drawn. In addition to the standard 
transmitted electron distribution, program JORDAN gives the distributions of both the 
small and large angle scattering. The program is given in Appendix D. With the 
“stand-alone" program several problem runs could be quickly conducted at many difl'er- 
ent energies. The results of those runs arc given in graphical form in Appendix E. l or 
the small angle distribution a discontinuity developed at energies less than 0.5 MeV’ at 
an angle of 105 degrees. At 0.01 MeV this discontinutity grows into a major peak. In 
this case the small angle distribution has been extended far beyond it's region of validity. 
We must remember that small angle approximatiotis were used in the dcvclopemcnt of 
the distribution equations. 

The large angle distribution' for the related energies show the same pattern. For 
energies less than 0.5 MeV the discontinutity at 105 degrees becomes large. One expla¬ 
nation is that the large angle scattering in the Jordan-Mack correction is based on a 
Possion distribution which requires the probability of a single event to be vciy small. 
As the energy decreases, the nutnber of large angle scatters becomes large and Poisson 
statistics arc not applicable. 
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VI. REMARKS 


This study looks only at the angular distribution of the transmitted electrons. Al¬ 
though the energy deposition for very thin foils is small, this aspect of the code and the 
correction should be investigated. Another of area of interest which needs elaboration 
is the apparent discontinuity in the small and large angle distributions at energies below 
0.5 .MeV. The problem occurs at 105 degrees and is not 2 dependent (See Appendix E). 
Finally, all the comparisons for this study were calculated using ITS version 2.0. An 
updated version 2.1 has recently been installed at the Naval Post Graduate School. Trial 
runs have shown the same problems in regard to very thin foils. Version 2.1 contains a 
correction to the Landau Straggling Algorithm and doesn't cflcct the subroutine AN¬ 
GLE. [Ref. I7j 
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VII. CONCLUSIONS 

The existing CLYTRAN code approximation for tlie multiple scattering in verv thin 
materials is inadequate when the thickness of the foil is less then three substep intervals. 
The Jordan-.Mack correction shows an improvement o\ cr a wide range of energies and 
materials. The Jordan-.Mack scheme does however break down for ver>- thin materials 
when the energj- of the incident electron beam is less than 0.5 McV. For moderate en¬ 
ergies the Jordan-.Mack method gives better results then the .Moliere first term approxi¬ 
mation. Of course if more MoUere terms were used this diflercnce would be minor. In 
addition, is was shown that although the large angle scattering distribution makes a 
small contribution, it should not be ignored for moderate energies. Future versions of 
the ns code system should include a correction to the multiple scattering for very thin 

materials. Finullv. the range of applicability must be determined and documented v.ith 
the final version. 
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APPENDIX A. ITS OV ERVIEW 


The Integrated Tiger Series of Electron Photon Transport codes (ITS) is the most widely 
used particle transport code in the world. The code package was developed to incorpo¬ 
rate eight individual codes wliich were developed over the period from 1968 to 1981. 
All the codes are based on the original ETRAN model developed by M. Berger and S. 
Seltzer. The ITS code system consists of four primary code packages (Ref. 1: p. 9): 

XD.AT.A: The electron and photon cross section data file. 

XGEN: The cross section generation program. 

ITS; I'he Monte Carlo program file. 

L PEML: .A machine portable update emulator. 

The heart ofthe ITS is the program library file ITS. which contains the eight Monte 
Carlo programs plus system directives for the CR.AY, IBM, VAX. and CDC operating 
systems. The update emulator program UPEML creates the various .Monte Carlo codes 
for a given system with any corrections to those codes that may be desired. The output 
fortran source code from UPEML is then compiled and stored as an executable module. 
Program XOEN generates the problem specific cross section data tape using file 
XD.AT.A for referenced inputs and a user defined input file. The Monte Carlo codes 
then read in the cross section tape and process the user defined problem. 

One ofthe eight ITS codes is CYLTR.AN, which simulates the transport of particle 
trajectories through a three-dimensional multimaterial cylinder. For this project only the 
CYLTRAN code was required. .As an ITS user the following steps were required to 
execute an ITS run 

1. Create the specific ITS code CYLTRAN with the required correction schemes. 

2. Generate a cross section tape based on the difi'erent t> pe of materials contained in 
the cylindrical geometry of a problem. 

3. Create an input file which list all the input parameters required to calculate desired 
outputs. 
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4. Submit the input file and the generated cross section tape to the ITS Monte Carlo 
codes to execute a run. 

Table 1 is a sample input file to generate the cross section tape for the materials in a 
cjlindrical geometr>-. Each material line represents a diffierent medium in the cylinder. 
Percentages of each material in a compound and its density must be specified. Single 
element lines such as Cu, has its density stored in ITS and is automatically used for the 
simulation when needed. 

Table 1. SAMPLE INPUT FILE TO CREATE A CROSS SECTION TAPE 


Energy H.fl 
Material Cu 
TITLE 

16 MeN’ Cross Section for Au foil 


Once a cross .section tape is generated, an input file with the parameters design Ibr 
a particular simulation must be created. Table 2 on page 41 is a sample input file to 
execute an ITS run. The ketxvords at the beginning of each line are relatively scli- 
explanatoiy. However, a few kexxvords and their numerical parameters require some 
explanation to fully understand their importance. These few kecAvords and tlieir expla¬ 
nation are listed [Ref. 1; pp. 2I-2?1: 

Geometry This kexAvord sets up the cylindrical geometry of the problem into input 
zones, where in this case there are seven input zones. Each line of numerical 
parameters following the keyword describes the dimensions of each zone. 

Electron/Photon-Escape This kexAvord tallies the number of incident electrons and 
photons that escapes the cylindrical geometry, either latterally or transmit¬ 
ted through. 

NBINE tallies the escaped electron photons in specified energy bins. 

NBINT tallies them in angular bins. 
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Electron/Photon-Flux This keyword tallies the energy deposition of electrons and 
photons in the subzones described in the kenvord geometry . 

Histories This kcyAvord represents the number of primary particle histories to follow 
as it transport through each medium in the cy linder. 


Table 2. SAMPLE OF AN INPinf FILE TO EXECUTE ITS 


Echo I 
TITLE 

Mat;Au 18.67 Mg/cm2 HistlOOOOO E* 15.77MeV 
Electrons 
Energy 79.0 
Cutoffs .1.1 
Position 0,(1 0,0 0.0 
Direction 0.0 0.0 
Geometrv I 

0.000 0.00096684 0.00 10.000 1 0 1 ] 

Electron-Escape 
NBINT 24 User 

0.5 I. 1.5 2. 2.5 3. 3.5 4. 4.5 6. 6.5 7. 

7.5 8. 8.5 9. 9.5 10. 15- 45.90. 180. 

NBINE 1 User 

0.1 

Histories 100000 


41 



APPENDIX B. JORDAN-MACK EQUATIONS 


The Jordan-.Mack correction to the multiple scattering for ver>' thin foils is based 
on a Williams Gaussian distribution. To determine the mean square angle, the screened 
Rutherford cross section with Moliere's screening angle are combined and integrated 
over the forward region: 


2 

O) 


re. 


Xt 


e^da 


•'0 


(5.1) 







P 


2( I - cos 
(1 — cos 0 + ;/)■ 


d{ cos 6) 


(5.2) 


where >/ = ' '.6\ is one half the value of Moliere's screening angle and the small angle 
approximation of O- = 2(1 — cos 6) has been used. The constants in this integral form 
the .Moliere unit probability angle, and with the substitutions n = cos - cos this 
equation siniplilies into a form which is easy to integrate: 
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This expression can further be simplified by the introduction of the constant /I 

2 _ >/ 

which is the form given in the theory section. To evaluate the 'argument" (I - + tj), 

the condition of one scatter occurring in the back region is imposed in the following 
form: 
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Rearranging yields; 
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APPENDIX C. FORTRAN SUBROUTINES 


A. ANGLE SUBROUTINE 

In the existing ITS code C^'LTRAN, the multiple scattering distribution near the 
escaping boundary' is calculated in a subroutine called ANGLE. As described in the 
theor}' section, this subroutine uses a Williams-type Gaussian approximation. In a 
subroutine called XINPLT, the GS average cosines with the corresponding D ranges are 
read from the cross section tape and stored in the variable names COS.AV and DRG. 
In the subroutine XPREP the variable name COSAV is transformed into a proportion¬ 
ality constant by the following code: 

DRGS(N,J) = DKG(N,J)/ISUB 

COSAV(N,J) = ( COKE-COSAV(N,J))/DRGS(N,J)XPREP 

where ISL'B is the number of substeps and CONE is the constant one. Since the D 
ranges are in units of Mg cm- this proportionality constant is then multiplied by the 
density in a subroutine called PREP so that the COS.AV values are ready for the AN¬ 
GLE subroutine. The crude approximation is done m only nine lines of code. When the 
proportionality terms COS.AV are multiplied by the distance to tlic edge boundary. 
SHD. the resulting term represents an approximation to the average cosine of the de¬ 
flection angle for that specific distance: 

ALF = COSAV(NT,MT)*SHD 
RA = RAN(IRA.N) 

IF (RA .LT. C1EH2) THEN 

ARG = MAX(-C88,-CTW0/ALF) 

EX = EXP(ARG) 

COM = CONE+ALF*LOG(RA*(CONE-EX)+EX) 

ELSE 

COM = CONE+ALF*LOG(RA) 

ENDIF 

The constants CONE, CTWO, C1EM2, and CSS are the values 1.0, 2.0, 0.002, and SS.O 
respectively. 


44 


Q0nf avQUciijit> to DTIC dc>«s not 

p0aBk hilly legible lepioductioo 


The expression lor the cosine of tlie dellected ande COM is dcrised from the M illiams 
Gaussian expression which is formed into a CDT and from which a random value is 
drawn: 
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B. ANGDET SUBROUTINE (JORDAN-MACK METHOD) 

For the Jordan-Mack correction scheme a subroutine called ANGDET was devel¬ 
oped to be called from the ITS subroutine ANGLE. All the previous coding just de¬ 
scribed for tlie crude appro.ximation was deleted by the following code: 


^DELETE ANGLE .XX-JX 


The call statement is then placed in the subroutine ANGLE with the patch subroutine 
•ANGDET placed at the end of ANGLE: 


*INSERT ANGLE.99 

CCN:iON/TE.Mr/JSUE ,I.MAT( INMT) ,!1PAIR,CTAX( 13) .NTAB,VT( INEM, LNIT), 
lAU'i I.NEM, EMIT) .2E( INEM,INMT) 

*I.NSERT A.NGLE. 120 

C. 

C AN’GLE DETERIJI NATION PATCH C 

C. 

CALL ANGDETt COM, SHD ,T.RH0( MT) , ZE( 1, MT;, ATVv’C 1, LT) ) 

C 

C.ENT OF PATCH.C 

’•'INSERT ANGLE. 190 

Sl'EROlTINE ANGDET( COM, SKD,E,RHO, 2, A) 

DOCBLE PRECISION COM,STH,CTH,SPH,CPH,ETA,FST,HST.AST, 

1 BST. >LMP, P1, THE2 ,TARG, SARG ,THE, ROT, PEXP, ALXMP, PCFN 
DOUBLE PRECISION R.RX.RF.RARG.ARG.COSMAIv.POF.POFZ.XTON 
DOUBLE PRECISION DSEED.IRAN 
COMMON /VAXRAN7 IRAN 
DIMENSION ROT{3.3),R(3,3),RX(3,3) 

DATA PI,R2ERC,EM2ER0,ALPHA,AVA/3. 1415926536D00,0. 281751 
1,0.511,137. 0371.0.6025/ 

RAN(DSEED)=GGUBFS(DSEED) 

Initialization of the Constants: 

ET*E/EMZERO + 1. 

B21» l./(ET*ET) 
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B2= 1.-B21 
23= Z**.333333 

ETA = .5*(B21/B2)*(Z3/(.885*ALPHA))**2 
1 *(1.13 + 3. 76*(Z/ALPHA)’M*'2/B2) 

FST = 2.*PI*RZERO**2 *Z*(Z+1) *B21/(B2*B2) 

1 *(AVA/A)*SHD*RHO 

Determination of the Integral Values: 

HST = 1. +ETA 
AST = 1. /(HST+1.0) 

RARG = 1./FST + l./(2.+ETA) 

ARG = 1./RARG 
COSMAX = HST - ARG 

COSMAX is the maximum angle such that one large angle scatter will 
occur in the back region. 

144 IF (COSMAX.LT.1.) GOTO 10 


When COSMAX is greater than one, a large angle scatter is not 
possible so that an explicit angle determination must be carried out. 

COSMAX =1.0 
BST = l./ETA -AST 
XMP = FST>''BST 
COM = 1. 

STH = 0. 

GOTO 1002 

When a large angle scatter is possible, the small angle portion 
is determined first using the Jordan-Mack small angle 
approximation equation: 

10 XMP = 1. 

BST * 1./FST 

PI = FST*(DL0G(ARG/ETA)-1. + ETA*RARG) 

THE2* 2 * PI 

TARG = 1. -DEXP(-PI*PI/THE2) 

The multiplication term SQRT(SIN(THETA)/THETA) as suggested by 
Bethe is samplied using a rejection technique: 

1000 SARG = 1. -TARG*RAN(IRAN) 

THE * DSQRT(-THE2*DLOG(SARG)) 

STH = DSIN(THE) 

REJECT = 1. 

IF (THE.GT. l.E-6) REJECT=DSQRT(STH/THE) 

IF (RAN(IRAN). GT. REJECT) GOTO 1000 

C0M=DC0S(THE) 

Initialization of arrays for the large angle portion; 

1002 R0T(1,1) = COM 
R0T(2,1) = 0. 




2030 

2040 

2050 


2060 

2070 


R0T(3,1) = -STH 
R0T(1,2) = 0. 

R0T(2.2) = 1. 

R0T(3,2) = 0. 

R0T(1,3) = STH 
ROT(2,3) = 0. 

R0T(3.3) = COM 

PFYP = -YMP 

ALXMP = DLOG(XflP) 

A random value is drawn from which the probability of each large 
angle scatter will be subtracted. If this random value is less 
than the the possion distributed value for one deflection, no 
large angle scattering will develope; 

POFN = RAN(IRA.M) 

XTON = 0. 

POFZ = DEXP(PEXP+XTON) 

IFCPOFN.LE.POFZ) GOTO 7777 


The random walk on the possion distribution is then carried out: 


POFN = POFN - POFZ 
DO 2000 N=l,500 
RF=RAN(IRAN) 

CPH = DMINl(COSMAX,DMAXl('l.D00,HST-l.D00/(AST+BST>>RF))) 
SPH = DSQRTCl.-CPH*CPH) 

THE = PI*(2.*RAN(IRAN) -1. ) 

CTH = DCOS(THE) 

STH = DSIN(THE) 

R(l,l) = CPH*CTH 
R(2,l) = CPH*STH 
R(3,l) = -SPH 
R(l,2) = -STH 
R(2,2) = CTH 
R(3,2) = 0. 

R(l,3) = SPH*CTH 
R(2,3) = SPK*STH 
R(3,3) = CPH 
DO 2050 1=1,3 
DO 2040 J=l,3 

RX(J,I) = 0. 

DO 2030 K=l,3 

RX(J,I)=RX(J,I)+R0T(J,K)*R(K,I) 

CONTINUE 
CONTINUE 
CONTINUE 
DO 2070 1=1,3 
DO 2060 J=l,3 

R0T(J,I) = RX(J,I) 

CONTINUE 
CONTINUE 
COM * R0T(3,3) 


The probability of another large angle scatter is then calculated; 
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XTON = XTON + ALXMP -ALOG(FLOAT(N)) 
POF = DEXP(PEXP+XTON) 

POFN = POFN - POF 
IF (POFN. LE.0.) GOTO 7777 

2000 CONTINUE 
7777 RETURN 
END 


C. ANGDET SUBROUTINE (MOLIERE METHOD) 

The Moliere approximation was programed into CYLTRAN in a similar mannor. 
The same delete and inserts for the subroutine ANGLE that were used for the Jordan- 
Mack method were used for the Moliere update. 


*INSERT ANGLE. 190 

SUBROUTINE ANGDET(COM,SHD,E,RHO,Z,A) 

DOUBLE PRECISION COM,STH,CTH,SPH,CPH,ETA,FST,HST,AST, 

1 B ST,XMP,P1,THE 2,TARG,SARG,THE,ROT,PEXP,ALXMP,POFN 
DOUBLE PRECISION R,RX,RF,RARG,ARG,COSMAX,POF,POFZ,XTON 
DOUBLE PRECISION DSEED,IRAN,C1,C2,BK,B,REJECT 
DOUBLE PRECISION CF,C3,A1,A2,A3 
DOUBLE PRECISION ALPHA,B2,Z 
COMMON /VAXRAN/ IRAN 

DATA PI,RZER0,EMZER0,ALPHA,AVA/3.1415926536D00,0.281751 
1,0.511,137.0371,0.6025/ 

RAN(DSEED)=GGUBFS(DSEED) 

Initialization of the constants: 

ET=E/EMZERO + 1. 

B21= 1. /(ET^'ET) 

B2= 1.-B21 
Z3= Z**.333333 

ETA = . 5*(B21/B2)*(Z3/(.885*ALPHA))**2 
1 *(1.13 + 3. 76*(Z/ALPHA)**2/B2) 

FST = 2. *PI*RZER0**2 *Z*(Z+1) *B21/(B2*B2) 

1 *(AVA/A)*SHD*RHO 


Evaluation of the integral values: 

HST = l.+ETA 
AST = 1./(HST+1. 0) 

RARG = 1./FST + l./(2.+ETA) 

ARG = 1./RARG 
COSMAX = HST - ARG 
BST = 1./FST 

Evaluation of Moliere constants: 

IF (Z.EQ, 13.D00) CF=-5. 2D00 
IF (Z.EQ,29.D00) CF=-5.6D00 
IF (Z.EQ. 79.DOO) CF=-6. 2D00 
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IF (Z.EQ. 82.D00) CF=-6. 3D00 
C2=6*FST/(7*ETA) 

C1=DL0G(C2) 

Al=l.ODOO 
A2=0.160D00 
A3=3, 33D00 

C3=(DL0G(A2*Z3*23*(A1+A3*((Z/ALPHA)**2/B2)))-CF)/(Z+A1) 

An alternate Jordan-Mack small angle approximation is done if 
the value of the constant Cl is less than one. This was done 
because Moliere's iterative equation diverges when Cl is less 
than 1.0. 

IF (Cl.GT. 1. ODOO) THEN 
BK=1.OIDOO 

444 B=BK 
BK=DL0G(B)+C1+C3 

IF (ABS(BK-B). GT.1.OD-4) GOTO 444 
THE2=2. 0D00*FST’’'B 
GOTO 446 

ENDIF 

445 THE2=2*FST*(DL0G(ARG/ETA)-1.ODOO+ETA*RARG) 

446 TARG = 1.ODOO -DEXP(-PI*PI/THE2) 

Again the multiplication factor is incorporated into the 
distribution: 

1000 SARG = l.ODOO -TARG*RANX IRAN) 

THE = DSQRT(-THE2*DL0G( SARG)) 

STH = DSIN(THE) 

REJECT = 1.ODOO 

IF (THE.GT.l.E-6) REJECT=DSQRT(STH/THE) 

IF (RAN(IRAN). GT. REJECT) GOTO 1000 
COM=DCOS(THE) 

RETURN 

END 
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APPENDIX D. FORTRAN SUBROUTINE JORDAN 


PROGRAM JORDAN 
DOUBLE PRECISION COM 
DOUBLE PRECISION IRAN.DSEED 
COMMON /VXRAN/ IRAN 
COMMON /DANS/ 

1 NDANl, NDAN2, NDAN3, NDAN4, NDANM 
DIMENSION DEG(100),CEG(100),HIS(99),FLX(99) 

RANNO(DSEED)=GGUBFS(DSEED) 

IRAN=63939 
NDAN1=0 
NDAN2=0 
NDAN3=0 
NDAN4=0 
NDANM=0 
NCNT=0 

OPEN (50,FILE='JINPUT',STATUS='OLD') 

OPEN (70,FILE='XYZ9') 

1 FORMAK/,/,/) 

WRITE (*,1) 

WRITE (*,*) ' 

WRITE (*,*) 

WRITE (*,*) ' 

WRITE (*,*) ' 

WRITE (*,*) ’ 

WRITE (*,*) ' 

WRITE (*,*) 

WRITE (*,*) ’ 

READ(*,*) NOPT 
READ (50,*) 

READ (50,*) 

READ (50,*) E,£ 

READ (50,*) 

READ (50,*) NW 
READ (50,*) 

RHOP=RHO*1000 
WRITE (70,*) 

WRITE (70,*) 

WRITE (70,*) 

WRITE (70,44) RHOP,Z,E,NUMHIS,NOPT 
44 FORMAT(F5.2,'G/CN2*,' MAT2:',F5.2,’ E:’,F5.2,'MEV HIS:',16, 

§ODEL:',11) 

NZERO=0 

NTWO=2 

NF0UR=4 

N0NE=1 

NTHR=3 

WRITE (70,*) 

WRITE (70,*) 

WRITE (70,*) 

WRITE (70,*) NONE,NONE,NZER0,NZER0 
WRITE (70,*) NONE,NZERO,NTVO,NZERO 


JORDANS PROGRAM' 

OPTIONS: 4. M0DEL4 (SMALL AND LARGE)' 

5. MODELS (MOLIERE)' 

6. MODELS (LARGE ONLY)' 

7. MODEL7 (SHALL ONLY)' 

ENTER: ’ 


D,RHO,Z,A 
IS,NUMBIN 
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WRITE (70,*) 

WRITE (70,*) 

WRITE (70,*) 

WRITE (70,*) NTHR,NONE,NONE,NONE 
WRITE (70,'(A8)') ' JORDAN' 

NUHBIP=NUMBIN-1 
WRITE (70,*) NUMBIP 
WRITE (70,*) 

IF (NUMBIN. LE. 0) GOTO 9999 

READ (50,*) (DEG(I),1=1,NUMBIN),DEG(NUMBIN+1) 

DO 1010 1=1,NUMBIN 

CEG(I)=C0S(DEG(I)*3.1415926536/180.) 

HIS(I)=0.0 
FLX(I)=0.0 
1010 CONTINUE 

DO 2000 N=1,NUMHIS 

IF (NOPT.EQ. 4) THEN 

CALL ANGDT4(C0M,SHD,E,RH0,Z,A) 

ELSE IF (NOPT.EQ. 5) THEN 

CALL ANGDT5(C0M,SHD,E,RH0,Z,A) 

ELSE IF (NOPT.EQ.6) THEN 

CALL AN'GDT6(COM,SHD,E,RHO,Z,A) 

ELSE IF (NOPT.EQ. 7) THEN 

CALL ANGDT7(C0M,SHD,E,RH0,Z,A) 

ENDIF 

DO 1020 1=1,NUMBIN 

IF (COM. GE.CEG(I+1)) GOTO 1030 
1020 CONTINUE 

I=NUMBIN 

1030 HIS(I)=HIS(I)+1.0 

IF(HIS(I).GT. HMAX) THEN 
HMAX=HIS(I) 

E.NDIF 

2000 CONTINUE 

XHIS=NUMHIS 
DO 2010 1=1,NUMBIN 

IF (HIS(I).LE.0.) GOTO 2010 

ERR=100.*SQRT((XHIS-HIS(I))/((XHIS-1.)*HIS(I))) 
DEGREE=(DEG(I)+DEG(I+l))/2. 0 
IF (I.EQ. 1) THEN 
DEGMIN=DEGREE 

ENDIF 

NER=ERR 

NCNT=NCNT+1 

WRITE (70,*) DEGREE,HIS(I),NER 
2010 CONTINUE 

WRITE (70,*) DEGMIN,DEGREE 
WRITE (70,*) NZERO.HMAX 

WRITE (70,*) ' VrVfjViVjVAiVALOOPS INFO**********' 

WRITE (70,*) NDANl 
WRITE (70,*) NDAN2 
WRITE (70,*) NDAN3 
WRITE (70,*) NDAN4 
WRITE (70,*) NDANM 
9999 STOP 
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END 

SUBROUTINE ANGDT4(COM,SHD,E,RHO,Z,A) 

DOUBLE PRECISION COM,STH,CTH,SPH,CPH,ETA,FST,HST,AST, 

1 BST,XMP,PI,THE2,TARG,SARG,THE,ROT,PEXP,ALXMP,POFN 
DOUBLE PRECISION R,RX,RF,RARG,ARG,COSMAX,POF,POFZ,XTON 
DOUBLE PRECISION DSEED,IRAN 
COMMON /VXRAN/ IRAN 
COMMON /DANS/ 

1 NDANl, NDAN2, NDAN3, NDAN4, NDANM 
DIMENSION ROT(3,3),R(3,3),RX(3,3) 

RANNOC DSEED)=GGUBFS(DSEED) 

DATA PI,RZERO.EMZERO,ALPHA,AVA/3. 1415926536D00,0. 281751 
1,0.511,137.0371,0.6025/ 

ET=E/EMZERO + 1. 

B21= 1. /(ET’VET) 

B2= 1.-B21 


Z3= Z**.333333 

ETA = .5*(B21/B2)*(Z3/(.885*ALPHA))**2 
1 *(1.13 + 3. 76*(Z/ALPHA)**2/B2) 

FST = 2. *PI*RZER0**2 *2*(Z+1) *B21/(B2*B2) 
1 *(AVA/A)*SHD*RHO 
HST = l.+ETA 


AST = 1. /(HST+1. 0) 

RARG = l./FST + 1./(2.+ETA) 
ARG = 1./RARG 
COSMAX = HST - ARG 
NDAN1=NDAN1+1 


144 IF (COSMAX. LT. 1. ) GOTO 10 
NDAN2=NDAN2+1 
COSMAX =1.0 


BST = 1. /ETA -AST 
XMP = FSr'BST 
COM = 1. 

STH = 0. 

GOTO 1002 
10 XMP = 1. 

NDAN3=NDAN3+1 • 


BST = 1./FST 

PI = FST*(DLOG(ARG/ETA)'l. + ETA*RARG) 
THE2 = 2.*P1 

TARG = 1. -DEXP(-PI*PI/THE2) 

1000 SARG = 1. -TARG*RANNO(IRAN) 

THE = DSQRT(-THE2*DLOG(SARG)) 

STH = DSIN(THE) 

REJECT = 1. 

IF (THE. GT. l.E-6) REJECT=DSQRT(STH/THE) 
IF (RANNO(IRAN).GT. REJECT) GOTO 1000 
CO.M=DCOS(THE) 

1002 ROT(l,l) = COM 
ROT(2,l) = 0. 

ROT(3,l) = -STH 
ROT(l,2) = 0. 

ROT(2,2) = 1. 

ROT(3,2) * 0. 

ROT(l,3) * STH 




5 


t 


!' 

R0T(2,3) = 0. 

■ R0T(3,3) = COM 

PEXP = -XMP 
ALXMP = DLOG(XMP) 

POFN = RANNO(IRAN) 

XTON = 0. 

POFZ = DEXP(PEXP+XTON) 

IF(POFN. LE. POFZ) GOTO 7777 
NDAN4=NDAN4+1 
POFN = POFN - POFZ 
DO 2000 N=l,500 

IF (N.GT. NDANM) THEN 
NDANM=N 

ENDIF 

RF=RANNO(IRAN) 

CPH = DMIN1(COSMAX,DMAX1(-1.DOO,HST-1.DOO/(AST+BST*RF))) 
SPH = DSQRTd. -CPH*CPH) 

THE = PI*(2.*RANN0(IRAN) -1.) 

CTH = DCOS(THE) 

STH = DSIN(THE) 

R(l,l) = CPH*CTH 

R(2,l) = CPH*STH 

R(3,l) = -SPH 

R(l,2) = -STH 

R(2,2) = CTH 

R(3,2) = 0. 

R(l,3) = SPH*CTH 

R(2,3) = SPH*STH 

R(3,3) = CPH 

DO 2050 1=1,3 
DO 2040 J=l,3 

RX(J,I) = 0. 

DO 2030 K=l,3 

RX(J,I)=RX(J,I)+ROT(J,K)*R(K,I) 

2030 CONTINUE 

2040 CONTINUE 

2050 CONTINUE 

DO 2070 1=1,3 
DO 2060 J=l,3 

ROT(J,I) = RX(J,I) 

2060 CONTINUE 

2070 CONTINUE 

COM = ROT(3,3) 

XTON = XTON + ALXMP -ALOG(FLOAT(N)) 

POF = DEXP(PEXP+XTON) 

POFN = POFN - POF . 

IF (POFN. LE. 0. ) GOTO 7777 
2000 CONTINUE 
7777 RETURN 

END wv VVHrWnWr* X.\AX.\X . *, 

SUBROUTI NT ANGDT5(COM,SHD,E,RHO,Z,A) 

DOUBLE PRECISION COM,STH,CTH,SPH,CPH,ETA,FST,HST,AST, 

1 BST,XMP,P1,THE2,TARG,SARG,THE,ROT,PEXP,ALXMP.POFN 
DOUBLE PRECISION R,RX,RF,RARG,ARG,COSMAX,POF.POFZ,XTON 
DOUBLE PRECISION DSEED,IRAN,C1,C2,BK,B,REJECT 
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DOUBLE PRECISION CF,C3,A1,A2,A3 
DOUBLE PRECISION ALPHA,B2,Z 
COMMON /VXRAN/ IRAN 
COMMON /DANS/ 

1 NDANl, NDAN2, NDAN3, NDAN4, NDANM 
DATA PI,RZERO,EMZERO,ALPHA,AVA/3. 1415926536D00,0. 281751 
1,0.511,137. 0371,0. 6025/ 

RANNOC DSEED)=GGUBFS(DSEED) 

ET=E/EMZERO + 1. 

B21= 1./(ET’'fET) 

B2= 1.-B21 
23= 2**.333333 

ETA = .5*(B21/B2)*(Z3/(.885*ALPHA))**2 
1 *(1.13 + 3. 76*(Z/ALPHA)**2/B2) 

FST = 2. *PI*RZERO**2 *Z*(Z+1) *B21/(B2*B2) 

1 *(AVA/A)*SHD*RHO 
HST = l.+ETA 
AST = 1. /(HST+1.0) 

RARG = 1./FST + l./(2.+ETA) 

ARG = 1./RARG 
COSMAX = HST - ARG 
NDAN1=NDAN1+1 
BST = 1./FST 
CF=-6. 2D00 
C2=6*FST/(7*ETA) 

Cl=DLOG(C2) 

Al=l.ODOO 
A2=0.160D00 
A3=3. 33D00 

C3=(DLOG(A2*Z3*Z3*(Al+A3*((Z/ALPHA)**2/B2)))-CF)/(Z+Al) 
IF (Cl.GT. 1. ODOO) THEN 
BK=1. OIDOO 

444 B=BK 
BK=DL0G(B)+C1+C3 

IF (ABS(BK-B).GT.l.OD-4) GOTO 444 
THE2=2. ODOO^-FSP^B 
GOTO 446 

END IF 

445 THE2=2*FST*(DLOG(ARG/ETA)-l. ODOO+ETA*RARG) 

NDAN2=NDAN2+1 

446 TARG = 1. ODOO -DEXP(-PI*PI/THE2) 

1000 SARG = l.ODOO -TARG*RANNO(IRAN) 

THE = DSQRT(-THE2*DL0G(SARG)) 

STH * DSIN(THE) 

REJECT * 1.ODOO 

IF (THE.GT. l.E-6) REJECT=DSQRT(STH/THE) 

IF (RANNOCIRAN). GT. REJECT) GOTO 1000 
COH=DCOS(THE) 

RETURN 

xx x xf!x fl xxxxxxxxx } . i i. xi . M .x } . A t.x Ax xx xxxx xx s ^ x i xnxxx xx x.\xxi . i x xx:.x 

SUBROUTINE ANGDT6(COM,SHD,E,RHO,Z,A) 

DOUBLE PRECISION COM,STH,CTH,SPH,CPH,ETA,FST,HST,AST, 

1 BST,XMP,P1,THE2,TARG,SARG,THE,ROT,PEXP,ALXMP,POFN 
DOUBLE PRECISION R,RX,RF,RAR6,ARG,C0SMAX,P0F,P0FZ,XT0N 
DOUBLE PRECISION IRAN,DSEED 


54 





COMMON /VXRAN/ IRAN 
COMMON /DANS/ 

1 NDANl, NDAN2, NDAN3. NDAN4, NDANM 
DIMENSION ROT(3,3),R(3,3),RX(3,3) 

DATA PI,RZERO,EM2£RO,ALPHA,AVA/3. 1415926536000,0.281751 
1,0.511,137.0371,0.6025/ 

RANNO(DSEED)=GGUBfS(DSEED) 

ET=E/EMZERO + 1. 

B21= l./(ET*ET) 

B2= 1.-B21 
23® Z**. 333333 

ETA * .5*(B21/B2)*(Z3/(.885*ALPHA))**2 
1 *(1.13 + 3.76*(Z/ALPHA)**2/B2) 

FST = 2. *PI*RZER0**2 *Z*(Z+1) *B21/(B2*B2) 

1 *(AVA/A)*SHD*RHO 
HST = 1. +ETA 
AST = 1. /(HST+1.0) 

RARG ® 1./FST + l./(2.+ETA) 

ARG * l./RARG 
COSMAX = HST - ARG 
NDAN1=NDAN1+1 

144 IF (COSMAX.LT.1. ) GOTO 10 
NDAN2=NDAN2+1 
CO.M ® 1. 0 
GOTO 7777 
10 XMP ® 1. 

NDAN3=NDAN3+1 
BST = l./FST 
COM=l. 0 
STH=0. 0 

1002 ROT(l,l) = COM 
BOT(2,l) ® 0. 

ROT(3,l) = -STH 
ROT(l,2) » 0. 

ROT(2,2) = 1. 

ROT(3,2) ® 0. 

R0T(1,3) * STH 
ROT(2,3) = 0. 

ROT(3,3) = COM 
PEXP ® -XMP 
ALXMP = DLOG(XMP) 

POFN * RANN0(IRAN) 

XTON = 0. 

POFZ = DEXP(PEXP+XTO.V) 

IF(POFN. LE. POFZ) GOTO 7777 
NDAN4«NDAN4+1 
POFN « POFN - POFZ 
NDANH«0 

DO 2000 N®1,500 

RF*RANNO(IRAN) 

CPH « DMIN1(COSMAX,DMAX1{-1.DOO,HST-1.DOO/(AST+BST*RF))) 
SPH » DSQRT(1.-CPH*CPH) 

THE » PI*(2.*RANNO(IRAN) -1.) 

CTH » DCOS(THE) 

STH a DSIN(THE) 

R(l,l) « CPH*CTH 
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2030 

2040 

2050 


R(2,l) = CPH*STH 
R(3,l) = -SPH 
R(1.2) = -STH 
R(2,2) = CTH 
R(3,2) = 0. 

R(l,3) = SPH*CTH 
R(2,3) = SPH*STH 
R(3,3) = CPH 
DO 2050 1=1,3 
DO 2040 J=l,3 

RX(J,I) = 0. 

DO 2030 K=1.3 


RX(J,I)=RX(J,I)+R0T(J,K)*RCK, 
CONTINUE 
CONTINUE 
CONTINUE 
DO 2070 1=1,3 


DO 2060 J=l,3 


ROT(J,I) = RX(J,I) 

2060 CONTINUE 

2070 CONTINUE 

COM = ROT(3,3) 

XTON = XTON + ALXMP -ALOG(FI.OAT(N)) 
POF = DEXP(PEXP+XTON) 

POFN = POFN - POF 
IF (POFN. LE. 0. ) GOTO 7777 
2000 CO.NTINUE 
7777 RETURN 


I) 


E.ND 


C***i>*iWfVr***,-t**i(HV***V,*,V*,WfVf*****,-r*,rtV*A,V*VnV**,V*****A,V,tVnV****>V**Vt*Vr*,WnV* 
SUBROUTINE ANGDT7(COM,SHD,E,RHO,Z,A) 

, ISION COM, STH,CTH, SPH, CPH, ETA, FST, HST, AST, 

PEXP, ALXMP, POFN 

PRECISION R,RK,RF,RARG,ARG,C0SMAX,POF,POFZ,XTON 
DOLBLE PRECISION IRAN.DSEED 
COMMON /VX'RAN/ IRAN 
COMMON /DANS/ 

1 NDANl, NDAN2, NDAN3, NDAN4, .NDANM 
DIMENSION R0T(3,3),R(3,3),RX(3,3) 

DATA PI,RZERO,EM2ERO,ALPHA,AVA/3.1415926536000,0.281751 
1,0.511,137.0371,0.6025/ 

RANNOC DSEED)=GGUBFS(DSEED) 

ET=E/EMZERO + 1. 

B21= l./(ET*ET) 

B2= 1.-B21 
Z3= Z*^». 333333 


ETA » .5*(B21/B2)*(Z3/(.885*ALPHA))**2 
1 *(1.13 + 3.76*(Z/ALPHA)**2/B2) 

FST = 2. *PI*RZERO’>*2 *Z*(Z+1) *B21/(B2*B21 

1 *(ava/a)*shd*rho 

HST ■ l.+ETA 
AST » 1./(HST+1.0) 

RARG « l./FST + 1./(2.+ETA) 

ARG * 1./RARG 
COSMAX » HST - ARG 
NDANl-NDANl+1 
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144 IF (COSMAX.lt. 1. ) GOTO 10 
NDAN2=NDAN2+1 
COSMAX =1.0 
BST = 1./ETA -AST 
XMP = FST’^^BST 
COM = 1. 

STH = 0. 

GOTO 7777 
10 XMP = 1. 

NDAN3=NDAN3+1 
BST = 1./FST 

PI = FST*(DL0G(ARG/ETA)-1. + ETA*RARG) 
THE2 = 2. *P1 


TARG = 1. -DEXP(-PI*PI/THE2) 
1000 SARG = 1. -TARG*RANNO(IRAN) 
THE = DSQRT(-THE2*DL0G(SARG)) 
STH = DSIN(THE) 

REJECT = 1. 


GT. l.E-6) REJECT=DSQRT(STH/THE) 
IF (RANNO(IRAN). GT.REJECT) GOTO 1000 
COM=DCOS(THE) 

7777 RETURN 


END 
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APPENDIX E. LARGE AND SMALL ANGLE SCATTERING 

DISTRIBUTIONS 


Since the Jordan-Mack correction scheme decouples the multiple scattering into 
large and small angle scattering, it seemed natural to plot the contribution of each sep¬ 
arately. The FORTR.AN program JORDAN was used to calculate the small and large 
angle scattering distributions. JORDAN is only a number tally calculation and involves 
only the columb scattering. Pair production, absorption, secondary’ production, and 
other particle interactions are not calculated. For veiy thin foils these efiects are minor. 
The input parameters to JORDAN were for the case of the Hanson gold foil, (18.67 mg 
cm-). The energy of the incident electrons was increased from 0.01 Mcv to 10.0 .\lev. 
The number of liistories was KK.1.1XMJ for all JORDAN runs. Figures IS thru 2-4 show 
the small angle distributions for this case. Figures 25 thru 34 show the large angle dis¬ 
tributions. Note the large discontinuity in the small angle distributions are at 1U5 de¬ 
grees. When the energy is greater than 0.5 MeV this discontinuity becomes scry nrinor. 








Figure IS. Small Angle Speclruni lor 0.01 .Me\: The iiipui parameters, were lor { 

a IS.66 mg cm- gold foil with 0.01 MeV intidem electrons. 1 



Angle (Deg) 






Small Angle Spectrum for O.U) Mc\: ilic mpm parameters, were the 
same as Figure 18 on page 59 except the incident electron enerev was 
0.10 .\IeV. 
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Figure 20. Small Angle ■Spectrum for 0.25 MeV: 1 lie input parameters were the 
same as Figure IS on page 59 except the incident electron energv was 
0.25 .\IeV. 
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Figure 21. Small Angle Spictruiii fur 0.35 MeV: Tlic input parameters were the 
same as Figure IS on page 59 except the incident electron energv was 
0.35 MeV. 
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Figure 2. 


1 . Small Angie Spectrum for 0.50 Me^■: The input parameters w ere ilic 
same as Figure IS on page 59 except the inciJeni electron enerev was 
0.50 McV. 
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igure 23. Small Angle Spectrum for 1.00 Me3’: 1 lie input parameters were the 
same as Figure IS on pace 59 except the incident electron enerev was 
1.00 McV. 
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Figure 24. Small Angle Spectrum lor 10.0 ,^Ie^'; The input par;inKiers woo the 
same as Figure IS on page 59 except the incident electron energy was 
10.0 McV. 
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Figure 26. Large Angle Spectrum for 0.02 MeV; Tlie input parameters were the 
same as Figure 25 on page 66 except the incident electron cnerev was 
0.02 McV. 
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Figure 27. Large Angle Spectrum for 0.05 ^le^’: The input parumciers were the 
same as Figure 25 on page 66 except the incident electron energy was 
0.05 MeV. 












Figure 28. Large Angle Spectrum for 0.07 MeV; The input parameters were the 
same as Figure 25 on page 66 except the incident electron energv was 
0.07 McV. 
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Figure 29. Large Angle Spectrum for U. 10 Me\': The input parameters were the 
same as Figure 25 on page 66 except the incident electron energv was 
0.10 .VleV. 
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Figure 30. Large Angle Spectrum lor 0.25 MeV: Tlic input parameters were tlie 
same as Fiaurc 25 on pase 66 except the incident electron cnerev was 
0.25 MeV.' 


71 












Figure 33. Large Angle Spectrum for l.OU Me\ : Tlie input parameters, were tlic 
same as Figure 25 on pace 66 except the incident electron energy was 
l.OU McV. 
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Figure 3*1. Large Angle Spectrum lor 10.0 MeN’: 7 lie input parameters were the 
same as Fisure 25 on page 66 except the incident electron energv was 
10.0 MeV.' 
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